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Abstract 

We find a closed-form for the distribution function (defined in terms of a 
Wigner operator) for hot cofoured particles in a background gluon field, in the 
hard thermal loop approximation. We verify that the current is the same as 
that derived from the known effective action. 



1 Introduction 

Thermal QCD is beset with IR divergence difficulties. A necessary first step 
is to sum the "hard thermal loops" — loop diagrams with external momenta 
small compared to the temperature T [1]. These hard thermal loops have been 
summed up to give an effective action [2]. In the hard thermal approximation, 
the hot particles (quarks or gluons) behave semi-classically; and this suggests the 
use of a transport equation. Such equations have been written down by treating 
the hot quanta as particles [3], and also derived by defining the distribution 
function as the thermal expectation value of a Wigner operator [4] . 

In the second approach, it has been usual to derive the transport equation 
from the Wigner operator, and then solve that equation in the hard thermal 
approximation. The purpose of this paper is to omit the intermediate step (the 
transport equation), and to find a closed form for the distribution function (that 
is, the expectation value of the Wigner operator) in the hard thermal approxi- 
mation. This closed form expression ought to encode the same information as 
the effective action. We verify this by showing that the current derived from the 
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distribution function is identical to that coming from the functional derivative 
of the effective action. 

In order to illustrate our method, wc first treat hot scalar "quarks" in the 
background colour field. Then we treat spinor quarks. Finally, we discuss the 
case where the hot particles are gluons, using the background field formalism. 



2 Hot scalar particles 

In order to illustrate our approach as simply as possible, we first treat the 

unphysical example of a hot scalar field (pa in the fundamental representation 
of SU{n), a = l,...,n, in a slowly varying background field A^^ab, where we 
take the n x n matrix A to be antihermitean and traceless and to include the 
coupling g in its definition. Wc work in Lorcntz space, with metric (+, — , — , — ). 

We define a distribution function to be the thermal expectation value (...) 
(with any zero-temperature part subtracted) of a Wigner operator. We use 
the symbol W to denote this expectation value. Thus we define (we denote 
hermitean conjugation by *) 
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(27r)4 
where 

and 



Wab{A;x,k) = 

d'^yexp{ik.y)Uac{A; x, x+){(t)c{x+; A)^*^{x-;A))[U*{A; x, X-)]db, (1) 

1 1 

x+=x+-y, x-=x--y, (2) 

.1/2 

U{A;x,x±) = Pe:iq)[ duy.A{x±uy)] (3) 
Jo 

and P denotes path-ordering of the matrices A. We have put A in the argument 
of in (1) to emphasize that (j) satisfies the field equation in the background 
field A. It follows from (1) that W is hermitean. 
Under a gauge transformation A, 

<l>{x) ^ A{x)-'<t>{x), (4) 

A^{x) ^A^ = A-'A^A + A-'d^A. (5) 



U{x, x+) A-\x)U{x, x+)A{x+y, (6) 



Then from (5) 
so that, from (1), 

W{x) A-\x)W{x)A{x). (7) 

To calculate W in (1), we need solutions of the field equation for (/>, in the 
background field A. This equation is 

[a^ + A^]2<^ = -mV (8) 
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(we give the scalar "quarks" a mass m, since this costs very httle extra comph- 
cation). We seek an approximate solution of the form 



X [K5(^; X, k) e^p{-ik.x)ab{k) + Vab{A; x, -k) eicp{ik.x)hl{k)] 



(9) 



where a is an annihilation operator and h* is a creation operator (normalized to 

[aa(k),a*(k')] = 2|k|<53(k-k')M- 

We assume that k is large compared to A, and to the rates of variation of V 
and A. Then (9) satisfies (8) provided that 



[2ik.d + 2ik.A - (5 + Af]V = 0. 



Write 



where (defining the covariant derivative D) 

{k.d + k.A)Vf) = k.DVo = 0, 

2ik.DVi « (a + AfVo = D^Vo. 
A unitary solution of (12) is 



Vo{A;k,x) = Pexp 



duk.A{x + uk) 



{k.d)-\k.A) + 



(10) 

(11) 

(12) 
(13) 

(14) 



Wc assume that k.d has a unique inverse when acting on functionals of A (this 
is discussed in equation (1.7) of [5]); so that we also have 



VQ{A;k,x) = Pexp 



duk.A{x + uk) 



= [Vo{A--k,x)]-\ (15) 



From (13) 



Vx » {2ik.D)-^D^Vo, 
« -V*D*^{2ik.D*)-^, 



(16) 
(17) 



where we make the convention that differential operators like D* act to the left. 

Using the solutions (9), and taking the thermal expectation value of the 
creation an annihilation operators, we get from (1) 
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Wab{x,k) = j dV'-^ j d^k'5+{k'''-m')N{\k'o\) 

[e-''-y[U{x,x+)V{x+,k')V*{x^,k')U*{x,x^)]ab + (fc' -fc')} , (18) 

where N is the Bose function, and the colours of the quarks in the heat bath 
have been summed over. 

Now lot B^{A\ , X, k) be defined to be the covariant field which coincides with 
A in the gauge where k.A = 0; so by definition 

k.B{A;x,k) =0. (19) 

This means that there is a gauge transformation, Vo{A; x, k) such that 

V^^B^Va = = V^^A^Va + V^^d^Va = V^^D^Vq. (20) 

Prom (19), this requires that Vq should satisfy (12), which justifies the use of 
the notation Vq. B satisfies the equation (using D defined in (26) below) 

k.VB^ = k.dB^ + k.AB^-B^k.A = k^F^^ = k.dA^-di^k.A + k.AAi,- A/^k-A. 

(21) 

This is right because it is a covariant equation, and in the gauge where k.A = 
it reduces to k.d{B — A) =0. 

In the hard thermal approximation, k is large compared compared to other 
dimensional quantities. We assume that this means that y in (1) is small com- 
pared to X. 

We will now separate three contributions to W. Type (i) is leading order, 
where we neglect y and Vi . This term is independent of A, and is simply 

W^i^ab = j^Ni\ko\)6{k-')Sab (22) 

as we would expect. 

In type (ii), we neglect Vi in (11) and work to first order in y in U and Vq. 
Type (iii) terms contain Vi, and for these we can neglect y in U and Vq. 

For the type (ii) term, then, we replace F by Vb in (18), and expand to first 
order in y, to get 

W(^^) = J^ J d'ye^'-y J d^k'6+{k'^-m'')N{\k'o\)[e-''''-yR{x,k',y)+{k' ^ -k')] 

(23) 

where 

R = ^y.A{x) + (^{y/2).^Vo^ V,* - Vo (^iy/2).^vi 

= y.B{x,k';A), (24) 
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where we have used (20) together with Vq^Vq = 1. 
So (23) gives 

= J^:^J^ik')N{\ko\)B^{x,k;A)], (25) 

where wc have used that B(—k) = B{k) from (21). 

In what follows, we will now find it useful to define a "2-sided" covariant 
derivative V, acting on matrices, by 

V^X = [D^,X] = -^X + [A^,X\. (26) 

If we let T) act on the product XY of two matrices (with the same argument 
a;), we get from (26) 

V^[X{x)Y{x)] = [D^X]Y + X[YDl], (27) 

where we make the convention that the conjugate operator D* acts to the left 
(as a differential operator), so that 

YV*^ = [D^Y*]*. (28) 

Note that, in spite of the notation, V is not an n x n matrix differential operator 
acting from the left. Rather it is an x matrix differential operator; so that, 
acting on a matrix function X, 

[DiiX]ab = {T^ij)ab,cdXcd- (29) 

Note that, from (21), 

k.V{B^, - A^) = -d^k.A. (30) 

We shall require the inverse of the operator k.V, which we will write shortly 
as {k.'D)~^. This is an x it? non-local operator. More explicitly: 

k.V{x,ab\{k.V)-^\x',a'h') = 6'^{x - x')5aa'Sbb' (31) 
and, for any matrix X{x), 

[{k.V)''XU{x)= j d^x'{x,ab\{k.V)-^\x',cd)X^d{x') (32) 

Now turn to type (iii) terms. In these, we may set j/ = in f/ and Vq. So, 
from (16), (17) and (18), we require S where 

2S{k') = [{k'.D)-^D''Vo\V^ - Vo[V^D*^{k'.D*)-% (33) 

where as usual D* is understood to act to the left. From (27), 

2k.VS = (r'Vo)yo* - Vq{V^D*'') = 2V^Y^ (34) 
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where 

2Y^ = {D^Vo)Vo* - Vo{Vo*D;) = 2B^, (35) 
(using (20)). Hence from (23) and using S{-k) = S{k), 

W^ui) = -^{k.V)-'V.B5{e)N{ko). (36) 

Combining the contribution from (36) with the type (ii) contribution (25), 
we finally get 



{2nf 



_d_ 

dk,. 



B^N{\ko\)6{e) ] + {k.V)-^V.BN6{e) 



(37) 



This is our result for the distribution function in the hard thermal loop approx- 
imation. 

The current is 

J^{x) =g j d^kk^W{x, k)=g j d^k0{ko)k^[W{x, k) - W{x, -k)] (38) 

being, like W, a,n n x n covariant matrix. Inserting (37) and integrating by 
parts, the temperature-dependent part of the current is 

Mx) = j^intd^kN{\ko\)S{k^)[B^ - k^{k.V)-^V.B]. (39) 

Note that 

= 0, (40) 

and J is hermitean. 

We now claim that this current is the functional derivative 

5V 

w- 



6At'{x) 

of an effective action (see equation (29) of [2]) 

r = j d'^x J d^kN{\ko\)6{k'MB{x, k; A)]\ (41) 

To verify this, we work out the variation of (37). From (21), this must have the 
form 

SB^{x)=dA^{x) + [f^k.dA]{x), (42) 
where / is some (nonlocal) functional of A. Therefore 

d[B^]=2[B^ + k^B.f]SA^. (43) 
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But since is an invariant, (43) must vanish for an infintesimal gauge trans- 
formation CO (using (5) and the definition (26)) 

6A^ - V^io; (44) 

so 

k.VB.f = -V.B (45) 

or 

B.f = -[k.V\-^V.B, (46) 

and thus (46) inserted into the variation of (41) yields the same as (39), as 
claimed. (Note that 6[iT:B'^] = -2Bab6Bab-) 

3 Dirac quarks. 

We now generalize from scalar quarks (j) to Dirac quarks tjj. Then (1) becomes 
W = J d^ye^^-yUii^m*, (47) 

where we suppress colour indices and spinor indices {W is a matrix in each sort 
of index). Instead of (9), we use 

i^aa = ^{27r)-''^'' I d''k6+(e - w?)[eM-^k.x)V::,^^4^ + ....], (48) 

where r labels two spin states, a, h are colour indices, a is a Dirac index. We 
expand 

KVa = y^abU^a^ + ^laM' (49) 

where 

(7.A; - m)u^''^ = 0. (50) 

Now act on ip in (48) with 
-(n.D + m){ij.D -m) = {j.D f =D^+m^ + ^i^^l^F^" (51) 

and use that {i'j.D — m)ip = 0. In a similar way to equation (13), we find that 
(suppressing colour indices and Dirac indices) 

2ik.DVl''^ = (D^Vo) + u^''^^{F^-'Voh^'y.u^''\ (52) 
When we work out W, there is a sum over the spin states r, and we use 

^uM^M = (A;.7) +m. (53) 



All the contributions, except the contribution from the F term in (52), will 
be like (37) except for an extra factor k.-j + m. 

In the contribution from the F term in (52), because the Vo and Vq are not 
differentiated, we may use VqVq* = 1. Thus we get (using (21)) 

2C=[jxlu{k.l + m)-{k.^ + m)jx7u]F^'' = -^F^^'k^lu = -2{k.V)B.^. (54) 
So finally, instead of (34), we have 
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[{k.V)-\V.B){k.^ + m) - B.j]N{\ko\)6{k^ - m^) 



— % 

:{k.'y+m) 



(2^) 

d 



-^{B^N{\ko\)5{e - m")] + {k.V)-^V.BN{\ko\)6{e - m") 



(55) 



(27r)3 

The current is now given by, in place of (38), 

Jfi,ab = 9 j d'^kivY^^Wab], (56) 

where the trace is over Dirac indices. Using 

tr[(fc.7 + TO)7^] =2fc^, (57) 

we get for J the same as (39) but with a factor 2 in front (there are 2 spin 
states). 

4 Hot gluons. 

We now turn to hot quantized gluons in a background classical Yang-Mills field. 
Call the background field A as in previous sections and the quantized field A, 
so that the total gluon field is A + A. Now we let the indices a,b,... stand 
for the adjoint representation. It may be convenient to think of Aab/^ as an 
(n^ — 1) X {n'^ — 1) antisymmetric matrix. We also define Ua by 

•^ab ~ dfabc^^c- 

The Lagrangian for a is 

^ — ^qu ~t~ ^nonqui (^^) 

where ^ 

C^qu = —-^Dabp.abv — DabvO-b)^ — a-a^iF^^I^ ttbv (59) 

is the part quadratic in a (which is the only part directly relevant to 1-loop 
digrams), and Cnonqu is the rest. Both parts of C are separately invariant under 

A^A-\DA), ^^A-MA, (60) 
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where as usual D = d + A. The invariance (60) wih guarantee that the W 
we derive transforms covariantly as in (6). On the other hand, the complete 
Lagrangian C is also invariant under 

A^K-^AK, K-^{D + A)K. (61) 

This gauge transformation justifies the fixing of the gauge in (58), for instance 
by adding the gauge-fixing term 

-\{D.af. (62) 

For the Wigner function, there are two possibilities (of the same structure 
as (1)): 

j (fy exp{ik.y){Ua^){x+){Ua,,){x-) (63) 



'"^ 2(27r) 

and 



^ / d*yexp{ik.y){Ufi^^){x+){Ufxp){x-) 



(64) 



were we have suppressed the group indices, and / is defined as 

//ii/ = ^^^lal, - V^a^. (65) 

(A Wigner function like (63) was intoroduced in [6].) In (64), m is some gluon 
mass, which is to tend to zero. 

Both (63) and (64) are invariant under (60), but (63) has the disadvantage 
that it breaks invariancx; under(61) ; so it may depend on the gauge-fixing of a. 
For (64), it is not clear how to carry out the limiting process consistently: it is 
well known that there is no smooth limit from a massive to a massless vector 
field alone (that is, without additional scalar particles). For this reason, we 
abandon W in (64) in this paper. 

In calculating the (expectation value of) the Wigner function, we need the 
spin sum of the polarization vectors e (analagous to (53)). In a "physical gauge" , 
we have 

E ^Ir^^^^Hfc) = -9,. + ^""^.^^^"^ = P,., (66) 

r=l,2 

where n is an arbitrary vector. If we could use W in (64), dependence on n 
would disappear. But using w in (63), a dependence on n may remain. 
The complete quadratic action (59) plus (62) may be written 

— {D^a^f -a^F'^'^a^. (67) 

The equation of motion is 

D^a^" - ^Fi^^a^ = 0. (68) 
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Like (49), we seek an approximate solution of (67) of the form 

Then, using (68) and steps like (12) and (13), 

i2k.VVl;^ = (2?Vo)eM + 2{F^,,Vo)e'^-^''. (70) 

The contributions from the first term (Vq) in (69) and the first term in (70) 
come out very similar to the scalar and Dirac cases, and give a contribution to 
w,. 



i 



2(27r)3 ^-^JBxN5{e)P^^] + ik.V)-'V.BNS{e)P^^'^ (71) 

where P^i, was defined in (66). The contribution from the "spin term", the 
second term in (70) is 

\{k.V)-^[nxFx,,k^ + kxFx,.nf,] + {n ^ ly) (72) 



2(27r)3 n.k' 

Neither (71) nor (72) can have physical significance in general, because of 
their dependence on the arbitrary vector n. However, the trace w = is 
independent of n: the trace of (72) is zero and the trace of (71) is similar to 
(37). We hope that we may use this trace with some confidence. So we take the 
current to be 

j^, = gj d^kk^w - g J d^k9{ko)k^[w{k) - w{-k)l (73) 

and then we get the expected contribution (39) from (71) (the gluon field is self- 
conjugate, but there are 2 polarization states: the trace of (66) gives = —2). 
In (73), we may us the property 

w{-k) = w'"'{k) (74) 

where cc denotes complex conjugation. Then all reference to negative values of 
fco in w is removed from (73). In order to prove (74), note that i appears in 
(63) and (10) only multiplied by k (since a and A are real). Since w* = w, we 
can alternatively (using (73)) write (74) as 

j^iab = 9 j d'^k9{ko)k^[wab{k) - Wba{k)]. (75) 



5 Conclusion 

The main results of this paper are the three closed expressions for the distri- 
bution functions: (37), (55) and (71) together with (72). In the last case, we 
expect only the trace to have physical significance. The secondary result is the 
verification that the currents obtained from these distribution functions coincide 
with the current coming from the known hard thermal effective action. 
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